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1. Introduction 

Recent developments in quantum technology have brought us a capability of 
manipulating and measuring a quantum system larger than two qubits using a number 
of different physical systems. The quantum states of such high-dimensional systems 
have been experimentally measured and characterized [1]. The standard method to 
evaluate quantum states in experiments is state tomography. With state tomography, 
we can reconstruct the density matrix of the system. Since the density matrix contains 
all the information of the quantum state we have, we are able in principal to calculate 
any characteristic of the system. The only problem is that the larger system gets, 
the exponentially more elements wc have to measure and compute to characterize the 
system, and the analysis of the quantum nature of a state quickly becomes intractable. 

Fortunately, there has been some recent progress in state characterization and 
visualization. A tomographic method for large systems with certain symmetries [2], and 
a state visuahzation method for discrete systems that extends the Wigner function [3, 4] 
have been developed. By contrast, for qunats (continuous variables) a long history exists 
for establishing a tool set that efficiently represents and analyzes quantum states in an 
infinite-dimensional Hilbcrt space. In particular, the standard method in quantum 
optics is based on mapping operator functions to corresponding c-number functions. 
For example, the Wigner function [5-9], Q-function [10], and P-function [11,12], which 
have been widely used in both theoretical and experimental analysis, are all c-number 
functions using coherent states. The coherent state is defined as {\a) = D{a,a*)\0) = 
^aa -Q;*a|g^-^ |--|^3j fjgj-g^ D{a, a*) is a displacement operator in the phase space where |0) 
is the vacuum. It can also be taken as the kernel that generates the Wigner function; 
hence [14] 

W(a, a*) = Tr[p • 2D(a, a*){-lf^D-\a, a*)], (1.1) 

Similarly, the Q-function and P-function can be obtained through (1.1) by replacing 
D{a,a*) with the appropriate operator, respectively. 

In this paper, we expand the use of these functions to general SU (N) systems by 
using 5" [/(A/") -symmetric coherent states. To begin, SU{2) coherent states for d-level 
systems can be generalized as the trajectory of the SU{2) group action U^^^ on the 
lowest wight state \ijjo) [16, 17]; hence 

|^,0) = C/f/^(^,0)|V'o) 

= E M + ^ ^^^id "'"'id e'(^-)^|-,m). (1.2) 

m=-f V 2 / 

Here an Euler angle decomposition was used for U^^'^{9,(f)) and d = M + 1, which is 
usually denoted as o? = 2j -|- 1 where j is a quantum number. We introduced the integer 
parameter M only for convenience in the generalization later in this paper. 

With the SU (2) coherent state given in (1.2), the Wigner function, Q-function, and 
P-function have been defined [18, 19]. They are successfully used to analyze atoms in a 



SU{N) -symmetric quasi-probability distribution functions 



3 



trap [19] and spin-squeezed states of two-component Bose-Einstein condensates [20, 21]. 
However, the main obstacle in using this state representation is the difficulty in adopting 
its composite structure into the analysis of the quantum system. For example, in 
quantum information processing, it is important to maintain the properties dependent 
on the composite structure of the system; quantities such as entanglement only have 
meaning with it. To accommodate a more detailed system structure, we first need to 
generahze the state representation method to SU (N) systems. 

Our starting point is to generahze (1.2) to SU{N). This generahzation can be done 



Here, = 6*1, 6*2, ... , O^^i, = 0i, 02, • • • , 4>n-i, and iV'o) is the lowest weight state in 
terms of the group operation Uj^ . The quantum number M defines the dimension d of 
the representation, hence the system size is given as 



We will revisit this and explicitly define the coherent states with an Euler angle 
parametrization [23] in Section 2. 

Using SU (3) coherent states, a Winger function has recently been constructed [24]. 
More general SU (iV)-symmetric distribution functions have been shown to exist [25, 26]. 
Furthermore, the Q-function can be generalized rather straightforwardly with these 
coherent state, however no general quasi-probability functions are constructed using the 
coherent states defined in (1.3). In this paper, we construct Wigner and P- functions for 
SU{N) in both the M = 1 and 2 cases. 

This paper is organized as follows. We first review the construction and properties 
of our generalized coherent states and then show how they help build a SU{N)- 
symmetric Wigner function, Q-function, and P-function through the Stratonovich-Weyl 
correspondence [15]. We then discuss some general properties of the functions and then 
conclude with an example. 

2. 5'C/(A'^)-symmetric coherent states 

We start by explicitly defining our SU{N) coherent states for rf-dimensional systems, 
where d is given in (1.4). First we denote the SU{N) generators [22] by the set 



as [22] 




(1.3) 




(1.4) 



{^N,M{k)} where A; = 1, 2, . . . , iV^ 
This set is made up of off-diagonal generators 



1. 



(2.1) 



A|^jj^(a, b) and A^^m(«' ^r a, b 



1,2,3,...,N; a <b, 



(2.2) 




(2.3) 
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A detailed procedure to construct these matrices and their properties is given 
in Appendix A. When c? = A^, i. e. M = 1, the representation is fundamental and the 
generators above reduce to the generalized Gell-Mann matrices {A^} for SU{N) [27, 28]. 
In particular, when d = 3, i. e. M = 1 and A" = 3, this procedure generates the standard 
Gell-Mann matrices {Xk} for SU{3): 



A«(l,3)^A3,i(4) 

A{J(2,3)^A3,i(6) 
and 

Ag>(3)^A3,i(3) = 



/ 1 \ 

1 , a£^(1,2)=A3,i(2) 



1 

I , Ag(l,3)=A3,i(5) 



1 

ODD 
1 
1 



,Ag>(2,3)^A3,i(7) 




(2.4) 



/ 1 \ 



v 



0-10 




,Ag>(8)^A3,i(8) 



1 

V3 



/ 1 \ 
1 
0-2 



(2.5) 



Given the generators {A7v,m(/^)}, we can employ the parametrization given in [23] 
for SU{N). Using this parametrization, a SU{N) operator for M can be decomposed 
as 



= ( n n AM{y,j{z))) X bm, 



where 



and 



N>z>2 2<y<z 



AM{y,j{z)) = e'^JV,M(3)</'(j,-i)+j(^)giAjv,M(l'2/)^(!/-i)+3(z)^ 



(2.6) 



(2.7) 



Bm ^ n e''^^^.*f((c+i)^-i)<^(jv(iv-i)/2)+c_ (-2.8) 

l<c<Af-l 

Here j{z) = for z = N and j{z) = ^^^^(A^ — i) for z ^ N . The Am (y, 7(2;)) terms 
are from the off-diagonal generators and the Bm term is from the diagonal generators. 
For example, for 5*^7(3), (2.6) gives us [29] 

u^u^{e,4>) (2.9) 
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Using (2.6), the generalized coherent state for SU{N) in the fundamental 
representation |(0,0)]y) can explicitly be written as 

(.ii<Pi+d2+-+4>N-2+<l>N-i) cos[6'i] cos[6'2] • 
_Qii-^i+<P2+-+<PN-2+<PN^i) sm[ei\ cos[6'2] 
_ei(03+<A4+-+0Ar-2+<Ajv-i) sin[6'2] cos[6'3] 



cos[6'Ar-2] sin[6'7v-i] 

■ ■ COs[6'Ar_2] sin[6'Ar_i] 

• cos[6'iv-2] sin[6'iv-i] 



+4'N-2+4'N-l) 



(2.10) 



sin[^^jv_4] cos[^^jv_3] cos[^^jv-2] sin[^jv-i] 
_gi(<^jv-2+0iv-i) sin[^^_3] cos[^jv-2] sin[^jv-i] 
_Qi4'N-i sin[6'Ar_2] sin[6'jv_i] 

COs[6'Ar_i] 

_. / 2(Af-l) , 

where cp — e y 'P(jv2+jv-2)/2 jg overall global phase. More generally, since we 
are looking at the lowest weight state as the reference state, SU {N) coherent states for 
d-dimensional systems can be easily shown to be equal to the ci-th column of [/j^(0, <^), 

\(e,ci>)^) = u^(0,ci>Mo) = [u^{e,<f>)U. (2.11) 

This is easy to see when {0, 4>) and |'0o) are represented in matrix form. The column 
vector of the lowest weight state has zero for all components apart from the bottom row 
element which is one. The only elements of the matrix U^f {6, cj)) that arc therefore 
relevant are those in the d-th cohimn. If, on the other hand, we had chosen the highest 
weight state as the reference state, the resulting coherent state would be the first column 
of U^{0,4>) i. e. [U^{6,4>)]i. Lastly, the coherent state (2.11) is equivalent to (1.2), 
regardless of the value of M, for N = 2 [16-18], as well as more general coherent 
states for larger values of N and M [22,30], if one makes the appropriate coordinate 
transformations on and </>. 

Lastly, by using the invariant volume element for the complex projective space in 
— 1 dimensions (dV^pjv-i) [31], 

= ( n ^(2/))d^^7V-id0^-i...dWi, 

2<y<N 

sin[2^i] y = 2, 

K(y) = { cos[^,_i]2^-3 sin[^^,_i] y^N, 
cos[^^_i]sin[^;v-i]''^-' y = N, 

derived from (2.6) we have the following resolution of unity for \{0, (f))\ 

2%'^ i(M!) Jg^^^_^ JeM-i,<l>N-i 
We denote the identity matrix of size d x d hy 1^ and we are using the following 
integration ranges [31], 



dVcpN- 



(2.12) 



(2.13) 



TT 



< < - and < 0i < 27r, for 1 < i < AT - 1, 



(2.14) 



to evaluate the integral. 
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As a final set of observations we note the following properties of our coherent states, 
(N + M-l) 



{{e,cpy^\AN,Mikm:(p)N)dVcpN-i = (2.15) 



27r^-i(M!) ' K-ut>.-i 

for all k, and, as a special case when M — 1, 

[ ■■■ [ ((0,0)jv|Aiv,i(a)|(0,0)]v) 



X ((0,0)]v|Aiv,i(6)|(0,</>)]v) dVcpN-i - 5ab, (2.16) 

where dab is the Kronecker's delta. 

3. 5't/(A^)-Symmetric Distribution Functions 

Having defined our generalized coherent states, we can immediately generalize the Q- 
function to SU (N) systems [20] . The Q- function of a density matrix p may be written 
down as 

Q(e,ci>)^{(e,cf>)^\p\(0,ci>)^). (3.1) 

Here, Tr[p] = 1 and, using (2.12), 
(A^ + M- 1)! 



zn [Ml) jg^^^^ JeN-i,4>N-i 



(3.2) 



Noting that 

{{6, <t>)^ \p\{e, ct>)^) = Tr[p . \{e, M){{e, 0)j^|], (3.3) 

we can see that we have the relation 

Q{e,cl>)^Tr:[p.\{e,<P)^){{d,<P)^\]. (3.4) 

Equation (3.4) is a special case of the more general Stratonovich-Weyl 
correspondence [15] that describes mappings between any Hilbert space operator X 
and a characteristic function /*(/3) on the classical phase space X {(3 = /3i, /32 ■ ■ ■ , Pn 
and 13 eX)hy 

n(3)^Tr[X.F%f3)]. (3.5) 

This mapping is very useful in that it allows us to represent a density matrix as a 
distribution function in phase space. The quasi-distribution functions are information 
complete to their original density matrix. This means that one can reconstruct the 
density matrix from its quasi-distribution function, /^(/3) h- )■ X, and the generating 
kernels F^{^) of /*(/3), which we will build from a set of hermitian generators (2.1), 
satisfy 

F^(/3) = F^(/3)tand / ... / F^(/3) di/(/3) = 1^. (3.6) 
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Here In is the iV-dimensional identity matrix and s determines the type of distribution 
function being described [32]: 

r^{(3) P-function, 
/+^(/3) ^ Q-function, 

f((3) ^ Wigner function. (3.7) 

Because of the way the /''(/3) are defined, they exhibit all the properties of the 
distribution functions they represent, however, from the application point of view, the 
value s should be chosen dependent on the properties to be investigated. For instance, 
the Q-function is easy to calculate, but often obscures the quantum nature in the states. 
In these cases, the Wigner function tends to be preferred to represent the signature of 
quantum properties by interference fringes. For finite-dimensional systems, (3.5) is a 
natural way to analyze the Wigner function, Q-function, and P-function, and we will 
build our generalized functions based on the various cases of (3.5). In detail, 

/wp(^'^) = ^[p-^^,m(^,<^)], (3.8) 

F^^}^{0,4>) is an operator of a M representation of SU{N), and (0, </>) denotes the 
parameters from the coherent states given in (1.3). 

Following (3.2) and (3.6), and using (2.12) as our integral measure, we will demand 
that fN,M,pi.^i4>) and F^^j^{e,cj)) satisfy 

(N + M-l)\ r f 

An [Ml) jfj^^^^ JeN_u<t>N-i 

and 

(N + M-lV f f 

/ITT [M\) Jg^^^^ JeN-l,<i>N-l 

For example, if we define F+\,(0,0) = \{e,cf))f}){{e,^)fj\ and /j^,\,,,(e, 0) = Q{d,(f>) 
then we can see that (3.9) and (3.10) are satisfied by (2.13) and (3.2). Lastly, to recover 
the density matrix, the relation 

^^ ^X^-W'?' / fl^:l{0,<P)Fii;M<l>)dVc,.-., (3.11) 

ZTT [Ml) Jq^^^^ J0N-i,<t>N-i 

has to be satisfied. 

As we have seen in the beginning, the generalized Q-function is rather easy to 
construct, however the Wigner and P-functions are somewhat more complicated. Thus, 
we will first look into the M — 1 case for arbitrary SU{N) systems. 

3.1. Fundamental Representation 

As before, we start with the Q-function, which, in the fundamental representation, can 
be written as 

/iU^,'^) = Tr[p-F+\(0,0)], (3.12) 
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where 

F+\(0,(/,) = |(0,(/,)jv)((0,0)jvl, 

= ^1^+ 1 E {i0,ci>yN\AN,iikm,^)N)^N,i{k). (3.13) 



k=l 

Now, in order for (3.12) to be useful, we need to be able to evaluate (3.11). We can 
only do this if we know F^\(0,0), which is the generating kernel for the P-function. 
Demanding that we satisfy (3.9) and (3.10), we see that F^\{e,(f)) IS 

1 + 1 

FN!iiO,<l>) = J^^N + ^ Yl ((^,0)]v|Aiv,i(A:)|(0,0)jv)Aiv,i(A;). (3.14) 

k=l 

With (3.14) done, the remaining function to generate is the Wigner function. Requiring 
that we again satisfy (3.9) and (3.10), we obtain for the Wigner function's generating 
kernel 

FNAe,<l>) = + Yl ((0,0)jv|Aiv,i(fc)l(^,0)]v)Aiv,i(fc)- (3.15) 

k=l 

Combining these results yields 

nAO,<t>)^-l^ + ^ Y ((0,</>)]v|Aiv,i(^)|(0,</>)]v)Aiv,i(^), 

k=l 

s = 0, 

n{s) = { 1 s = +1, (3.16) 

s = -1. 

Substitution of (3.16) into (3.8), with M = 1, yields results that agree with existing 
Wigner, Q-, and P-functions for iV = 2 and iV = 3 [24,33]. Lastly, (3.16) with (3.8) 
satisfies (3.9), (3.10), and (3.11) for all three s values. 

Next, in the fundamental representation, a density matrix p can be represented 
by [34] 

1 /iV^ 

P = ]^liv + V E ^k^NAk). (3.17) 

k=l 

For N — 2, the set of all pure states are characterized by n • n, n = {rife}, while for 
N > 2, the set of all pure states are characterized by n ■ n and n t*: n = n with the star 
product being defined in (A. 2). By substituting (3.17) into (3.8) and using (3.16), as 
well as exploiting (A. 3), we get 

fhA^,<t>) = ^ + ]^^j^m E ((^,0)]v|A^,i(^)l(^,'^)]vK. (3.18) 

k=l 
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To recover the density matrix from this function, we substitute (3.16) and (3.18) 
into (3.11), as well as exploit (2.15) and (2.16), to get (for the s = case) 

m 

p 



2tt^ 
1 



T / ■ ■ ■ / <l>)K,iiO, <t>) dVcp.-., 



N 



2N 



For the pure state case, the density matrix pps 
decomposition [35], 



2{N-i) ^^l^^-^(^)l^)- 



(3.19) 

has the following n 
(3.20) 



Substitution of (3.20) into (3.18) and (3.19) yields 
and 



J2 {{e,4>)i,\AN,i{km,<p)N)mN,i{km, 



k=l 



Af2-1 



PPS = j^^N + l Yl (V'|Aiv,i(A:)|V')Aiv,i(A:), 

k=l 



(3.21) 



(3.22) 



which is equivalent to (3.13). Similar calculations can be done for the s = ±1 cases. 
3.2. Higher Dimensional Representations 

This problem cannot be simultaneously resolved for all three s values by simply replacing 
\{6, 0)]v) with more general spin coherent state representations from (1.3), or like those 
defined in Refs. [22,30], the full form of the kernels must be calculated. To accomplish 
this, we start by generahzing (3.16) to generate a M = 2 representation of F^^j^{d, 0), 

1 ^ 

n,2{0:(l>) = -,td + J2uj%^2{c>N,2{^,(l>,c), (3.23) 



c=l 



where, using (1.4) and Section 2, 

(6Ar,2-(c-i))^-l 

U^,2{d,<l>,c) = J2 ((^''^)^IV2-(c-i),c(/^)|(0,0)^)V,_(,_,^,c(/c), 



(3.24) 



k=l 



and 



1 



'N+2,2 



1 



bN+2, 



^N,2(c) = < 



2&Ar+2,l 
1 

2 

c = 2, 



2u 



^,2(1)) 



c=l, 
c = 2. 



1, 



s = 0, 



s = +1, 



(3.25) 
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It can be shown that (3.23) with (3.8) satisfies (3.9), (3.10), and (3.11) for all three s 
values. 

As an example, we note that for a M = 2 SU (2) system, when s = 0, our new 
FN,2i.^^(t>) gives us 

^2,2(^,0) = ^Id + $^C^2%(C)^^2,2(0,0,C), 



c=l 



1 * 

= -I3 + J2 <l>)l\A3,i{km, <t>)l)A,,i{k) 



(3.26) 



fe=i 



+ a;°,(2) J] ((0, </,)^|A2,2(fc)|(0, (/>)^)A2,2(A;), 

fc=i 

. g . — . 2 

= + 4^ ((0, </>)^|Afe|(0, ct>)l)\k + ^"^^"3"^^°^ ((0, <t>)l\Jk\{e, ct>)l)Jk. 

k=l k=l 

Here we have again used (1.4) and Section 2 to recognize that A2^2{K) = Jk, the spin-1 
representation of the SU{2) Pauli spin matrices, and that As^i{k) = Xk, the standard 
Gell-Mann matrices for SU{3). Lastly, evaluating (1.3) for N,M — 2 yields 



-2i4>2 



e^'-^i sin[^i] 

siri[26'i] 



e-2i^i cos[^i]2 

It is easy to verify that (3.26) yields an equivalent operator for the Wigner function as 
that given in Refs. [33, 36, 37] and, using (3.8), gives equivalent Wigner functions as that 
given in [19]. Furthermore, a similar calculation for the s = ±1 cases yield equivalent 
operators as that given in Rcf. [33] for the Q-function and P-function. Therefore, despite 
the apparent differences between expressions, (3.23) is the appropriate generalization 
of (3.16) to M = 2. 



V2 



(3.27) 



4. Discussion and Conclusion 

In this section, we mainly discuss two issues: the graphical representation of states of N- 
level systems using our distribution functions and the correspondence relation between 
the various distribution functions. First, we show some examples of the graphical 
representation for a SU (4) system by considering a Werner state [38] of two qubits, 

/ 1-7 \ 
_ 1 1 + 7 -27 

PWerner T „ o i i n 

4 —27 1 + 7 

y 1-7 / 

Here the parameter 7 defines the purity of the state. In particular 7 = 1 corresponds a 
pure state whereas 7 = gives the completely mixed state. Evaluating pwemer via (3.17) 



(4.1) 
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gives us 




(4.2) 



which can be substituted into (3.18) to yield 



/4,i,Pwe_(^' = ^{^ - + 4cosM - (1 - 6cos[2^i] cos[^2]' 



- 3cos[2^2] - 12cos[(/)i - 02] sinful] sin[2e2]) sinf^g]^ )Q(s) 



(4.3) 



Here, n{s) is as defined in (3.16) with AT = 4, i. e. 0(0) = \/5, ^^(+1) = 1, and 
fi(-l) = 5. 



The function given in (4.3) is expressed in a five-parameter space, thus it is not easy 
to represent its entire property in one figure. So, we take cross sections of the function. 
To do this efficiently, we look at the element of the phases 0i and 02 in (4-3), that is 
12cos[(;/)i — 02]- This element is efi^ectively one parameter {(pi — <^2), so we can set it to 
have two extreme cases: 0i, 02 = and 0i = tt, 02 = 0. Doing this gives 



Figures 1 and 2 show the parameter regions where these quasi-distribution functions 
exhibit negative values for both cases 0i, 02 = and 0i = tt, 02 = respectively. As we 
expect, the Q-function docs not yield negative values in the entire parameter regime, 
however when the states are pure enough, the P-function and the Wigner function can be 
negative. When the purity of the Werner state, indicated by 7, becomes small enough, i. 
e. when the state is more mixed, the parameter regime for the negative value disappears. 
In fact, the P-function can be negative when p> 1/4, while the Wigner function shows 
a negative region when p > 1/2. Such negativity in these quasi-distribution functions is 
often considered evidence of quantum nature in the states, however, our results indicates 
such a simple explanation does not apply to SU (N) systems. 

There are of course other ways to reduce the dimensionality of the parameter space. 
The optimality of a representation is dependent on the properties we are after. Taking 
a cross section is the easiest way to generate a graphical representation to discern rough 
properties of the state. 

Next, we briefiy discuss the correspondence between different distribution functions. 
As we mentioned before, the concrete expression of the various distribution functions 
depends on the parametrization of the SU{N) group operators. In this paper, we 
employed the parametrization given in [23], which is an extension of the ones used 
for SU{2) and SUiZ). This parametrization gives us a way to write out the Wigner 
function, Q-function, and P-function that, in the M = 1 and M = 2 cases, makes 
their correspondence easy to see. Furthermore, it shows us how more general M 
representations should be related. 



2 + 4cos[2^3] - (1 - 6cos[2^i] cos[^2]^ 




(4.4) 
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Q function 




Figure 1. Graphical representation of the Q (top), P (middle), and Wigner (bottom) 
functions from (4.4) for 4>ii4'2 = and various values of 7. The three plots at far left 
show the case of the completely mixed state with 7 = and the far right ones plot 
the pure state with 7 = 1. The value of 7 for each plot is at right of the graphic. The 
colored areas represent those values of 9i, 62, and 9^ where the corresponding function 
is negative. 



To begin, for the M = 1 case, we start by making the following redefinition of (3.16), 
F^,,(0,0) = F^,i(0,0)-^ll^, (4.5) 



such that 



FUe,cj.) = ^F^^,{e,cj.) (4.6) 



is true. This can only be done because our Q{s) terms in (3.16) are all positive definite. 
Using (4.5) and (4.6) we therefore get 

f'^,i(e,*) = ^f';,i(«» + ^iA'. 

From this we can see how we can convert between the various ^{6, (p) operators and, 
through (3.8), the various fN,i,pi^y 4^)- 



SU{N) -symmetric quasi-probability distribution functions 13 



■Q function 




Figure 2. Graphical representation of the Q (top), P (middle), and Wigner (bottom) 
functions from (4.4) for </>! = tt, 02 = and various values of 7. The three plots at far 
left show the case of the completely mixed state with 7 = and the far right ones plot 
the pure state with 7 = 1. The value of 7 for each plot is at right of the graphic. The 
colored areas represent those values of 9i, 62, and 9^ where the corresponding function 
is negative. 



A similar procedure can be done for the M = 2 case. In detail, following (4.5) we 
redefine (3.23) to give us 

F^,2(^,0) = n,2(^,0)- Jld- (4.8) 

By construction, ^^2(^5*^) is non-singular and of rank d. It is therefore invertible, 
allowing us to have 

F^;2(e, 0) = 0) ■ FN,2i^, ■ FnA^, 4>)- (4.9) 

This allows us to state (via (4.8)) that 

F^;^,{e, 0) = F^;^,{e, ct>) ■ F^,2(0, 0)-i (f^,2(^> 0) - ^i"^) + \^'^^ 

= T(s', s, 2) ■ F^,2(^, 0) + (1 - T(s', 2)) ■ ^1,, (4.10) 

where we have made the following definition: T(s', s,2) = F^ 2(^)0) ' F^ 2i^^'P)^^- 
Conversion between the various f'^2pi^' 4^) '^i^ (3.8) is now straightforward. 

In general, we can see that the transformation sequence between the various 
functions in the general M case is equivalent to (4.10) if we make the following definition: 
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T(s', s, M) = F^'m(^, </>) • F^,m{^, <t>)~^- For example, when M = 1, T(s', s, 1) reduces 
to Q{s')/Q{s) as expected. 

To conclude, in this paper we have given an explicit set of S'f/(A^)-symmetric 
functions that represent finite-dimensional versions of the Wigner function, Q-function 
and P-function by using generalized coherent state. In the case of the general M 
SU{2) and M — 1 SU{3) representations, these functions are equivalent to previously 
derived finite-dimensional Wigner, Q-, and P-functions with an appropriate parameter 
change [24, 33, 36, 37, 39, 40]. The quasi-probability distribution functions in this paper 
have been generalized to a higher quantum number M — 2. Such quasi-probability 
distribution functions may also have some benefits to characterize qubit-qunat systems. 
These hybrid systems are becoming extensively investigated in the context of quantum 
information devices. For more complex systems, there are possibilities to further 
generalize the formula to an arbitrary M. However, the analysis has showed that the 
process is not as straightforward as the SU (2) case [19] and further work will be necessary 
to complete the generalization. 



Acknowledgements 

We would like to thank Jon Dowling, Bill Munro, and Peter Turner for helpful 
discussions. This work was partly supported by NICT and JSPS. 



Appendix A. 

The {Ajv,m(^)} matrices given in (2.1) are a subset of the corresponding Lie algebra 
of SU{N); a set of Hermitian, traceless matrices of size d x d that are defined in the 
following way [22] : 

(i) Define a general basis |mi, m2, . . . , mAr) where M — X^feLi^^fe; "^fe ^ ^) 
M e Z+. 

(ii) Define the following three operators: 

Jh\mi, . . . ,ma,mb, . . . , thn) = \/ {ma + l)m;, |mi, . . . , -M, m;, - 1, . . . , mjv) 
for 1 < a < 6 < A^, 

J^\mx, . . .,mb,ma, ■ ■ .,mN) = \^ma{mb + 1) |mi, . . . ,ma - l,mb + 1, . . . .tun) 
ioY 1 <h < a < N, and 



J^|mi, . . . , mc, . . . , mAr) = a / ^ (S^ irik - crric+i] |mi, . . . , mc, . . . , tun) 

yc(c+l)V^ ) 

for 1 < c < - 1. 

(iii) Using the basis given in (i) and the operators given in (ii), define the following 



matrices: 

,{1} h\ ^ , Tb 



A}^Ua,b)^J^ + Jl 
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Ag>^(a,6)^-i(J,"-J„^), 

ASUc+ir-l)^J^ (A.l) 

for a, 6 = 1, 2, 3, . . . , iV; a < 6 and c = 1, 2, . . . , iV - 1. 

(iv) Combine the three matrices given in (A.l) to yield the set {Ajv,m(^)} where 
= l,2,...,iV2-l. 

In general, our lambda matrices can be used to define the M = 1 star product, 



i^*y)k = y ^1^,2)2 Tr[{A^,i(0,A^,i(j)}-A^,i(fc)]a;,|/,-, (A.2) 
as well as be shown to satisfy 

[^N,M{i),-^N,Mij)] = C X fijk^NM^^ 

fijk = ^ X T>[[Aiv,M(i), Aiv,M(j)] • ^N,M{k)], (A.3) 

thus forming a basis for the corresponding vector space, and a representation of the spin 
generators of SU{N). For example, when M — 1 and c = 2i, (A.l) and (A.3) reproduce 
the form, and properties, of the generahzed Gell-Mann matrices [27, 28] . 
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